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1 Introduction and overview 



In this paper we investigate the role that dissipation can play in biasing va- 
cuum selection after a symmetry breaking phase transition. Our work was 
initially motivated by a mechanism suggested by Nagasawa and Branden- 
berger P to stabilise non-topological classical solutions by out-of-equilibrium 
effects. They studied an 0(4) model in which charged and neutral scalars 
coupled differently to a thermalised photon bath. Under the assumption that 
only the charged fields receive thermal corrections from the bath, these au- 
thors argued that the vacuum manifold would collapse from to 5*^ and 
therefore temporarily stabilise non-topological strings. 

The analysis in [1 is focused on the stability of embedded strings when 
immersed in a thermalised plasma. However, the primordial question of 
whether the formation of the defect is dynamically favoured has never been 
addressed. By looking at the requirements that favour their formation we 
found a close link between "external" dissipation and vacuum selection that 
ranges beyond our initial aim. By external we mean a source of dissipation 
that is not related to the interactions between the system and the heat bath. 

The vacuum selection we discuss can take place in the early universe and 
in heavy-ion collisions. In the early universe, the most relevant areas for 
applications are in the studies of preheating at the end of inflation |2j and 
defect formation in non-equilibrium cosmological phase transitions 

Current and future heavy-ion collision experiments also provide scena- 
rios where dissipative vacuum selection might take place. A specific problem 
where this process might occur is the formation of disoriented chiral conden- 
sates jUini- Processes of similar nature might take place in the quark-gluon 
plasma where modes with different thermalisation times coexist. However, 
our study, which is based on a symmetry breaking, does not provide by itself 
a mechanism for this case. 

The mechanism we investigate is illustrated here for the same 0(4) scalar 
field theory in 3 + 1 dimensions used in Ij. Some general assumptions are 
required to specify the properties of the model but when presented out of a 
particular context, at first sight, these might not appear natural. Therefore, 
we think it convenient to start with a brief outline of our study. In the context 
of our work we see the field theory as an effective model for the soft long- 
wavelength modes of a system coupled to a heat bath. At low temperatures 
the system has a symmetry broken phase and the symmetry is restored above 
some finite temperature T^. 
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The following are key requirements concerning the coupling of the scalar 
fields to the heat bath. In the ordered phase, two of the scalar fields, say 
for definiteness 0i^2, have decoupled from the heat bath, while the remaining 
fields, say ^3,4, stay coupled. This situation is more natural than it seems at 
first sight. For example it occurs when the fields 0i,2 are neutral and ^3,4 are 
charged, both with respect to the same conserved charge, and the heat bath 
consists predominantly of quanta of the associated gauge field. 

Furthermore, the coupling between ^3^4 and the bath is assumed to be 
much stronger than the scalar self-coupling. Then, we expect that by the 
time 03,4 thermalise at the heat bath temperature T < Tc, the 0i,2 fields are 
still out of equilibrium. In other words, the relative strength of the couplings 
implies that the relaxation times of ^3,4 are much smaller than the ones of 
01,2- This is a condition that can be realised in the early universe. The less 
favourable transition is probably the most recent one, the chiral symmetry 
breaking transition, because of the strong pion interactions. Under these 
circumstances, when some fields are coupled and other uncoupled to a heat 
bath, the question then arises: What is the effective vacuum manifold in this 
system? 

When all the fields are coupled to the heat bath and dissipate according 
to the fiuctuation-dissipation relation the answer is well known. The vacuum 
manifold is the three-sphere covered by all the equivalent scalar configu- 
rations that minimise the free energy. The soft modes described by the fields 
do not form a closed system because their energy is being exchanged with the 
heat bath, e.g. via collisions and decay channels, giving origin to dissipation. 
Let us assume, as it is normally the case, that these terms are characterised 
by viscosity coefficients rji associated with each scalar 0j. 

The situation we study here differs in two ways. First, the neutral fields 
are not coupled directly to a heat bath. Second, and most importantly, 
we consider these decoupled fields to have "external" sources of dissipation. 
For example, in a cosmological context this type of dissipation is naturally 
associated with the expansion of the universe. As a result the decoupled 
01,2 fields stabilise in steady states that can be characterised by an effective 
temperature T^s < T, independently on whether these fields condense or not. 

The most interesting effect occurs when the "external" dissipation is much 
smaller than the dissipation in the coupled fields due to their interaction with 
the bath. We would expect a small external dissipation to have only a negli- 
gible effect on the evolution of the fields. However, this is in general not the 
case. We have that at vanishing scalar self-coupling the limit 771,2 — > 0, i.e. 
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the "external" sources of dissipation are "switched off" , is singular. This has 
the effect of changing the symmetry breaking pattern even for small values 
of ?7i,2- In this limit the neutral fields are selected to condense, therefore ef- 
fectively reducing the vacuum manifold from to as originally suggested 
by Nagasawa and Brandenberger [I]. 

In reality, the uncoupled fields are still receiving energy from the bath 
and the role of the external dissipation is simply to release this energy at a 
comparable rate. This rate is much smaller than the corresponding values 
for the coupled fields and has no appreciable effect on them. 

The origin of the vacuum manifold reduction in our analysis is identified 
without having to call for out-of-equilibrium effects as in We emphasise 
the important role played by the existence of different steady states for the 
various fields due to the "external" source of dissipation. 

The vacuum selection takes place above a small critical dissipation which 
occurs when the neutral fields stabilise at a "cold" enough Tes- Our conclu- 
sions, therefore, is that the reduction is more widely applicable than previ- 
ously thought. 

Our simulations are governed by phenomenological Langevin equations 
describing the dynamical evolution of the fields. These equations have been 
previously used in a relativistic context to study non-equilibrium phenomena 
in cosmological phase transition [HI 13 13 • There are known limitations to the 
use of these equations which we discuss in Sect. 4 and 5. They provide 
nevertheless an economic and qualitative good description of the different 
processes involved in the dynamic evolution of the fields where the coupling 
to the heat bath is expressed by rapidly fluctuating fields and the dissipation 
effects are expressed by viscosity terms. In particular, this makes it easy to 
analyse the effects of dissipation terms that are not related to interactions 
with the heat bath. 

The paper is organised as follows. In Sect. 2, we discuss the role of 
decoupling and dissipation in the vacuum biasing mechanism. A toy model 
illustrating the dissipation requirement in a solvable system is presented in 
Sect. 3. We illustrate in this simple case the effects of non thermal dissipa- 
tion terms that are not related to exchanges with a heat bath but originate 
for example from the expansion of the system. In Sect. 4 we study dissipa- 
tive vacuum selection in the 0(4) symmetric model in 3+1 dimensions. This 
section is divided in four subsections. We begin with a discussion on the use 
of phenomenological Langevin equations. The results of our simulations are 
then presented in Sect. 4.2 and in Sect. 4.3 we analyse the effect of varying 
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the parameters of the model. Finally, in Sect. 4.4 we discuss the condi- 
tions that guarantee the vacuum selection to be in place when all fluctuating 
(and dissipative) contributions are accounted for. We end in Sect. 5 with a 
summary and a discussion of future work. 

2 Decoupling and dissipation 

We analyse the field dynamics in a system that undergoes a symmetry break- 
ing transition and where different fields sectors in the system reach distinct 
"thermalised" states with some heat bath. The nature of the bath will be 
characterised below. Some fields arrive at a standard thermalised state at 
the temperature of the bath after a relatively short relaxation time. The 
remaining fields stay out of equilibrium for a longer period, which can still 
be small compared to observation times. The fields that take longer to either 
thermalise or reach another type of equilibrium state are weakly or indirectly 
coupled to a heat bath. We refer to them as decoupled in a convenient loose 
sense that will be made clear later. In particular, we are interested in the 
situation where the decoupled fields condense following a finite temperature 
phase transition. In order for this to happen these fields must lose most of 
their energy and for this reason we will follow closely the role of dissipation. 

In order to discuss a setting where this scenario can be realised we use 
an 0(4) linear sigma model Lagrangean, 

C^^d,cf>,d^cf>i -^(cf>icP,-v'Y, ^ = 1,2,3,4, (1) 

to describe the propagation and self-interactions of soft modes. For conve- 
nience, we consider (pi and 02 to be neutral scalars and 03 and 04 to be the 
constituents of a charged scalar 0=*" = (03 ± i(\)/^j\pl with respect to a f/(l) 
charge. For concreteness let it be U (1)em- Therefore, 0^ is coupled to a bath 
of photons while the neutral scalars are not. In a more comprehensive ana- 
lysis the effects of the fiuctuations from the hard modes of the scalar fields 
are also to be taken into account. We will discuss their effects in the next 
sections. 

By keeping the model simple we can aim at a better understanding on 
how a small dissipation can play a role in selecting the vacuum. This is the 
effect we wish to emphasise and alongside we lay the conditions under which 
the vacuum manifold shows the selectiveness that favours the formation of 
such embedded structures. 



4 



Throughout this paper we work under the assumption that the couphng 
between the charged scalars and the photon bath is much stronger than the 
scalar self-couphng. Without this assumption leading to an effective sepa- 
ration of scales no non-trivial selection seems to take place. With these 
conditions, we expect the charged scalars to thermalise "quickly". Their 
relaxation time sets the scale for what we will refer as a quick thermalisa- 
tion time. In practice, at the observation scale, the charged scalars can be 
said to remain in equilibrium. The neutral scalars have of course a longer 
thermalisation time and are at least close to thermal equilibrium. 

One of our aims is to express quantitatively the distinction between the 
steady state reached by the neutral and the charged fields. If the decoupled 
fields have no direct process to dissipate their energy they reach a thermal 
equilibrium state at the temperature T of the photon bath. This therma- 
lisation occurs because the decoupled fields are not completely cut off from 
the photon bath due to the quartic scalar self-interaction. The rate at which 
the neutral fields thermalise depends on the strength of A. Elsewise, if they 
dissipate due to the expansion of the system as it cools, as in the early 
universe or heavy-ion collisions, the steady state they reach is "colder" . This 
effect leads to their selective condensation and suggests the use of an effective 
temperature Tes as a way to parametrise the distinct steady states. We will 
present a detailed discussion of T^s in the next sections. 

Before ending this section we refer the reader to the further underlying 
view that the Lagrangean (Q) is better suited for the broken phase. Our 
model can be seen as a sector in a larger theory. For instance above Tc 
other mediating bosons that decouple at the transition can be responsible 
for maintaining all the scalars in thermal equilibrium. A known example 
that inspires this view can be found in the electroweak transition. In this 
case the W- and Z-bosons become massive below the transition while the 
photons remain massless through the transition. 

Renormalisation due to radiative and thermal corrections coming from 
hard modes could also have been included in (jT)) but this should not have 
a major effect in our phenomenological analysis. More important are the 
fluctuation-dissipation effects coming from the scalar hard modes. They 
do indeed have some effect in the main simple picture laid down in this 
section but they do not change the outcome if the relatively magnitudes of 
the couplings remains as assumed above as we show at the end of the next 
section. 

Before we address the dynamical description of the model described by 
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(0) we look at a simple two particle system where we can perform analytic 
calculations that highlight the main features of the "vacuum biasing" . 



3 Biasing in a two field system in zero dimen- 
sions 

In this section we illustrate how decoupling and dissipation combine to bias 
the distribution of energy in a toy system. We have two coupled one-dimen- 
sional oscillators of unit mass and restoring constants ttIq and m^. Alter- 
natively, they can be interpreted as two interacting fields in zero spatial 
dimensions. The system can be seen as a first approximation for an effec- 
tive model for the zero Fourier components of a field theory in three spatial 
dimensions. It is a simplified version of the 0(4) model introduced in the 
previous section, but now with only two fields. Here one "charged" field 0c 
is coupled to a heat bath, and the other "neutral" field 0o is not coupled. 

We adopt a Langevin description, in analogy to what we will do in the 
next section for the system described by (P), where the coupling to the bath 
is represented by a random rapidly varying field. In a canonical form the 
equations of motion are, 

00 = TTo 

TTo = -rio-Ko - mQ(pQ - — 

O(po 

4>c = rcc 

dV 

The interaction between the fields is contained in the potential V{(f)Q, 0c) and 
the dissipation is expressed by the viscosity terms involving the coefficients 
?7o and rjc. The interaction of the charged field, (pc with the "photon" bath 
is modeled by the random field ^, 

m) = o, mm) = mt-t'), (3) 

where Q is the variance of the Gaussian white noise. This ensemble of fields 
can be described by the probability density p(0o, ttq, 0c, vTc, t). The evolution 
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of this distribution is governed by the Fokker-Planck type equation, 



dp dp f 2 dV\ dp d{nop) 

+ (ml^. + ^)^ + + . (4) 



;/ diTc duc 2 dri 



For any physical observable A(0o, ttq, 0c, t^c) we have 

{A){t) = d(l)odTTod(l)cd7rc A{(l)o,7ro,(j)c,'n:c)p{4>o,'rco,(j)c,7rc,t), (5) 



and its time derivative can be calculated using equation (jU and integrating 
by parts. The equation of motion for (A) then reads 

d(A) , dA^ ,dV dA^ 2,, dA^ , dA^ , dA^ 
dV dA^ ^ dA^ , dA^ n ,d^A^ 

In equilibrium we have dt{A) = for all physical observables A. One way of 
obtaining a full description of the system in equilibrium is to obtain all the 
expectation values of all combinations of powers of the four quantities (pQ, 
ttq, (pc and Tic- These are the A^-point functions for (classical)-fields in equili- 
brium. It is easy to see that except for the simple case of quadratic potentials 
(where we expect the system to be solvable), these equations mix correlations 
of powers of different degrees in the dynamical variables. This leads to an 
infinite number of linked equations similar to the Dyson- Schwinger hierar- 
chy. Here we will consider only the quadratic case and take the interaction 
potential to be 

V^(0O,0c) = A0o0c. (7) 

This choice will enable us to carry an analytic treatment. For the case with 
a quartic interaction X(pl4>l we have checked that the main features of the 
results are the same. 

For the potential (|7j), the second order correlation functions with i ^ j 
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satisfy the following system of 10 equations, 



71^ \{-Ki(t)j) + m2(0i7ri) + ri,{-Kf) - 6^2^/2 = 

0i7rj -> A(0,^) + m]{(j)i(j)j) + T]j{(f)i'Kj) - {-KiTij) = , (8) 

which is closed because it involves only expectation values of degree two 
in the dynamical variables. After some work we can solve it obtaining the 
average for the momentum of the decoupled field, 



(^0 + n 

In the simpler case when mo = rric = m this reduces to 



cJ 



The other momentum is related to this one by 

(.?)^|;-|(.„^). (10) 

Next we assume that the dissipation for the coupled field and the amplitude 
of the noise obey the fluctuation-dissipation relation, Q = 2r]cT, where T 
is the temperature of the bath. We now define the effective stabilisation 
temperatures 



n^ = {^!) = T-^K^). (11) 

Vc 

for the neutral and charged field respectively. In the case of zero dissipation 
coefficient for the neutral field, r/o = 0, not only these temperatures are the 
same but also T^q = T^'g = T for any finite non vanishing values of A. On the 
other hand, when r/o > the effective temperatures of the asymptotic states 
of each field are different which justifies using the quantities defined by (jllj) 
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as convenient parameters to express quantitatively the distinct steady states. 
We note that similar definitions have been introduced in glassy systems jH], 
however the Tes we use here are constants as they refer to steady states after 
relaxation instead of long lasting transient states. 

In this toy model the biasing mechanism corresponds to an unequal dis- 
tribution of the kinetic energy between the different types of fields. The main 
features that characterise it are easy to identify from equations (jllj) . To start 
with, note that if rjo <^ rjc, T^'^ is hardly affected with relation to its value at 
?7o = 0, whereas T^g can be noticeably reduced, provided < m^rjoric. When 
these two conditions are simultaneously satisfied even a very small value of 
the ratio rj = 77o/'7c can cause a large effect on the asymptotic configuration 
of the uncharged fields 0o- We have verified that simulations for a model 
with quartic self-interactions indicate the existence of analogous relations. 

The reason why a small value of rj is able to lead to a qualitatively different 
behaviour is the existence of a singularity to which we now turn our attention. 
Consider the two sequences of limits, A ^ followed by r^o ^ and its 
reverse when 7]q is taken to zero first. In the former we have Tgg — * and 
Tgg — *• T, and in the latter both T^g and T^'^ approach T. Clearly, the 
phenomena that we are studying corresponds to a regime along the first 
sequence of limits where the two types of fields are weakly coupled to each 
other but a dissipation in the neutral field prevails even after the self-coupling 
is "switched-off" . Therefore, we are not only considering t] small but also A 
has to be "small" in the sense that 

<V = ^«l. (12) 
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We might expect that the physical origin of the sizeable effect for small rj 
follows from r/o not being associated with fluctuation-dissipation effects. But 
this by itself is not sufficient to explain what we observe. Equally important 
is the existence of at least two scales in the problem. In order to clarify this 
point let us consider the case when the neutral fleld also evolves under the 
effect of fluctuations. To this end we modify the second equation in (j21) which 
now reads 

^0 = -^0^0 - - jr- + Co, (13) 

O(po 

where the dissipation coefficient rjo is now decomposed into two terms, t]q = 
^ext _|_ ^fl rpj^g term corresponds to an external source of dissipation 
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while the second one is related to fluctuating forces according to a standard 
fluctuation-dissipation relation, 

ieoit)) = 0, ieom^it')) = ^o^it - f), (u) 

where Qq = 2riQTQ, with Tq the temperature of the bath coupled to the neutral 
fleld. It is natural in a flrst approximation to take Tq = T as we should 
not expect a selective behaviour for the high energy modes. Moreover, this 
situation corresponds to a case where biasing is less favourable. The strength 
of the coupling of each fleld to the heat bath reflects itself in the relative 
magnitudes of the dissipation coefficients. This matter will be made clearer 
in the full 0(4) model. 

A generalisation of the calculation leading to including the noise for 
the neutral fields and with Tq = T gives 



T' —T 

eff eff '/o 



ext 



-/(mVo/A'), (15) 

70 

with f{x) = x/(l + x), from which we see that the effective temperatures 
of the steady states can still differ substantially if in addition to (fT^ we 
have ?7o^* ~ 7]^. These combine to guarantee that the argument of / is large 
therefore assuring that / ~ 1. Under these conditions T^'g — T^g ~ T. 

Note that an external dissipation should for consistency affect both fields. 
The inclusion of such non-thermal dissipation in the equations of the charged 
fields adds a new term in the right-hand side of (fT^ , analogous to the present 
one, where for instance the prefactor is now rj^^^/rjc instead of rj^^^/rjo- As 
^ext _ ^ext ^j^jg iga^^g small, O{rjo/rjc), reduction in the thermalisation 
temperature of the charged fields and no qualitative change. 

Therefore, it is justifiable, under the requirements at the end of the last 
paragraph together with (jisp . to neglect the external dissipation contribution 
to the charged field as r]"^^^ ~ rj^^^ -C r]^. It is the combination between 
the presence of an external source of dissipation and the existence of two 
scales, say t]^ and r]^, that is behind the physical origin of the effect we are 
analysing. The effect is only noticeable when the external dissipation is at 
least comparable to the smaller scale. In the next section we carry out a 
more explicit discussion on the role of these scales. 

We conclude the study of this toy model with a couple of graphical illus- 
trations of the distinct steady states and their associated relaxation times. 
For convenience we go back to the starting study case, = and r/g = '7o''* 
with the system governed by equations 0. 
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Figure 1: The effective stabilisation temperatures for each of the fields as a 
function of the viscosity coefficient rjQ. The curves shown are for three values 
of A and m = T = rjc = 1. 

Equations are represented graphically in Figure 1 for three different 
values of A. Clearly, as rjQ increases from zero, T^s decreases from the starting 
value T. The drop in temperature is more pronounced when the self coupling 
is smaller. Teg approaches negligible values for rjo > rjc and sufficiently small 
A < 0.2. This simply reflects that when the indirect interaction of 0o with 
the heat bath is weaker this field dissipates its energy more efficiently. On 
the other hand, the field 0c because of its direct coupling to the heat bath 
remains "hot", i.e. Tg^ stays close to T, but the larger A is, the more it 
deviates from T. This deviation is a natural consequence of the direct dissi- 
pative viscosity term for 0o which is not balanced by any fluctuating effects. 
Although not shown in Figure 1 we have from (fTT|) that when A — oo the 
effective temperatures become identical for any value of tjq. 

The relaxation times for the decoupled fields for various values of the self- 
coupling are shown in Figure 2. The results were obtained by determining 
the eigenvalues of the homogeneous version of equations (j21) {i.e. ignoring 
the noise term). The several relaxation time scales for the system are pro- 
portional to the inverse of the real part of these eigenvalues. In Figure 2 we 
show for each value of rjo the largest of these time scales, which we interpret 
as the equilibration time for the uncoupled fields. We also checked that for 
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^0 

Figure 2: The relaxation times for (pQ as a function of the viscosity coefficient 
rjo. The curves are for three values of A and m = T = rjc = 1. 

all values of rjo and A there is one eigenvalue leading to a relaxation time 
close to 2.0, corresponding to the charged fields. When the two values of 
the dissipation are the same, for rjo = rjc = I, the two time scales coincide 
as expected. The rapid slowing down in the rate at which 0o stabihses as A 
decreases shows a dependence closer to A~^ as we would expect. This naive 
expectation derives from the effective noise term 0{X) that the interacting 
potential induces in the equation for 0o in the equations of motion 

In summary, by looking at a two field toy model we observe an unevenness 
in the way the kinetic energy is distributed between the "charged" and the 
"neutral" field. The conditions for this biasing are possible because of the 
presence of an external dissipation term and different strengths for the cou- 
pling of the fields to the heat bath. When the field with the weakest coupling 
to the bath is decoupled we recognise that this counterintuitive behaviour is 
due to a singularity in the limit of vanishing external dissipation. In a realis- 
tic situation the weakest coupling to the bath should not be neglected. The 
singular gives place to a two scale regime and the biasing is expect to occur 
when r/^^* extracts energy at least at a rate comparable with the input from 
the bath coming from the weaker coupling. In the next section we will see 
how a similar effect contributes to a vacuum selection following a symmetry 
breaking phase transition. 
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4 Vacuum biasing in the 0(4) model 



4.1 The use of the Langevin approach 

We now return to the model described by the Lagrangean (PJ) and we will 
study its dynamical evolution under the Langevin approach already used in 
Sect. 3. In order to simulate the dynamics of the fields we use phenomeno- 
logical Langevin equations 



with = Af ^, V being the T = vacuum expectation value, and where rji 
and are respectively the viscosity coefficients and the Gaussian noises. For 
the fields that couple to the photon bath and that are assumed to thermalise 
at its temperature T, we have 



where (3VLi = 2rii, with /3 = 1/T, according to the fiuctuation-dissipation the- 
orem. Below Tc this relation is used only for i = 3,4, while in the disordered 
phase it is assumed for all the fields. 

In the ordered phase we consider the neutral fields to be decoupled from 
the bath, = 0, but let r]i ^ 0. The non vanishing value of these coefficients 
are due to an "external" source of dissipation. An expansion of the system, 
as in the early universe, is a possible origin for this type of dissipative term. 
In fact, this has been the main motivation behind this study. The toy model 
in the previous section showed us that even a small external dissipation could 
give rise to non negligible effects. In cosmology we expect a similar situation 
(except at very early times) in the sense that the dissipation due to the 
expansion of the universe is much smaller than the one associated to the 
charged ffelds which is due to the interactions with the "photon" bath. 

It is in place to say something on the advantages and limitations of using 
the Langevin equations (fT^ at this stage. These equations describe the 
classical out-of-equilibrium evolution of a system of coupled ffelds. Therefore, 
at most, it provides effective equations for the long- wavelength modes of 
quantum ffelds with large occupation number. For this reason it is often 
used as a phenomenological set of equations to study close-to-equilibrium 
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(16) 



k=l 



{Ux,t)^j{x',t')) = nAjS^^\x - x')6{t - t'). 



(17) 
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effects near pfiase transitions ^U]. Tfie use of tlie relativistic form of these 
type of equations has been motivated by cosmological apphcations [HI Ej . 

The relative simphcity of these equations is their main asset. We take 
the noise to be white and the dissipative kernel local both in space and 
time. Nevertheless, they provide a good starting set of equations to look for 
qualitative answers to questions in close-to-equilibrium dynamics provided 
one keeps a careful perspective of its shortcomings. This is the approach we 
take here. 

The shortcomings of equations ()16|) are best understood when a more 
systematic derivation of the effective out-of-equilibrium dynamics from first 
principles is carried out . The phenomenological approach to the 

study of hot non-Abelian plasmas provides also useful insights jElEl- There 
are, for instance, the questions on whether the origin of the noise is external 
or internal p!3], or that the Langevin approach is only reliable near thermal 
equilibrium when a quasi-particle type of approximation can be justified. Let 
us say something briefly on the second issue. 

One of the safest ways of understanding the regimes where a Langevin 
description applies in quantum field theory is to start from the Kadanoff- 
Baym equations |16j and work out under which conditions these equations 
can be interpreted as the result of Langevin processes |T7]. The high tem- 
perature limit is a well established condition and is naturally implicit in (fTBj) 
as this is an effective equation for long-wavelength modes, k <^ T. The use 
of classical equations is also justifiable by the more recent investigations on 
the reliability of the classical field theory limit to the dynamics of quantum 
fields out of equilibrium [THl dl El at high T or near Tg. 

4.2 The simulations 

We use a discretised version of ()16|) in three dimensional square lattices with 
50'^ to simulate the evolution of the 0(4) model (P). A leap-frog algorithm 
with time step 6t = 0.05 is used. Larger lattices of 100^ have been used to 
verify the stability of our results. A Gaussian random number generator is 
used for the rapidly changing fluctuations. All quantities are measured in 
units of the T = vacuum expectation v. The dimensionless quantities are 
identified with a tilde. For example, v = 1 and f5 = v/T. When choosing 
the lattice spacing we need to ensure that the modes with wavelength longer 
than ~ are not cut off. In our runnings merely for reference we took 
physical scales from the chiral symmetry breaking effective mean field model. 
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By using a lattice spacing Ax = 0.25 for = 93 MeV we can work up to 

temperatures of approximately T ~ 350 MeV. 

We run our simulations for successive temperatures of the heat bath de- 
termined by r2+ = 2?7+T, the amplitude for the noise of the charged scalars. 
Starting from a temperature T > Tc we bring down the temperature across 
Tf.. While in the disordered phase all the fields are taken to be in contact 
with the heat bath. It is only for T < Tc, when the vacuum expectation 
value starts to increase, that the neutral fields are decoupled. 

In Figure 3, we plot the order parameters for the condensates of neutral 
and charged scalars against the inverse temperature p. These are, respec- 
tively, (101021) and (103041), where (|0i0j|) = \jY.n=i,jiv Iv are ave- 
rages over the entire lattice. As the details of the decouphng of the neutral 
fields are not known we show the curves for three different types of decou- 
pling. The label of instantaneous, exponential and power law refer to the 
way the variance Qq = ^lo{T) approaches zero starting from its value at the 
disordered phase as T continues to decrease below Tc. For the runnings in 
Figure 3, the values of all the viscosity coefficients rji are kept the same for 
all the scalars. The different decouplings do not have an effect on the final 
values of the order parameters which always favour a non vanishing value for 
(|0i02|). This shows a bias for the condensation of the neutral fields resulting 
from the relative large value of the ratio between the viscosity coefficient of 
the neutral by the charged scalar. Here we used rjo/r)^ — 1 to emphasise the 
case when the condensation of the neutral sector is strongly favoured. 

For the large value of rjQ used in the simulations for Figure 3 the neu- 
tral and the charged scalars have an equally effective channel to cool as the 
temperature decreases. However, the fiuctuations coming from the interac- 
tion with the heat bath have the effect of slowing down the dissipation of the 
charged fields. This will favour the neutral fields to roll down more effectively 
to the bottom of the potential and condense. 

The neutral fields are not blind to the photon bath due to the scalar self- 
coupling. They dissipate through the viscosity term but they gain energy 
via the scalar self-couphng. However, as long as the effects of fluctuations 
hitting the neutral modes is small and any increment of energy can quickly 
be dissipated, which occurs when A is not too large and tjq not too small 
compared to rj^, the neutral scalars continue to monopolise the vacuum. As 
we discuss next the situation might change as rjo decreases. 

The curves in Figure 4, where we set 770 = but all other parameters are 
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Figure 3: Diagram showing the variation in the order parameter as the system 
is cooled across T^.. We use A = 1 and a 50^ box. We show curves for three 
different decouphngs: instantaneous, exponential and power law. 

kept the same as in the runnings used for Figure 3, show that in this case the 
charged scalars condense. This suggests that as t^o is decreased below some 
critical value the charged fields condense instead of the neutral ones. This 
results in a "superconducting" background where the photons are massive 
which clearly violates our working assumption of having a thermalised photon 
bath. Therefore, we conclude that the condensation of the charged scalars 
is an artifact of our simulations in this region of the parameter space of the 
model. 

From the simulations we have discussed, we anticipate the existence of 
a critical value for r] = tiq/ti^ in the interval < rjcr < I- A precise 
determination of r/cr is numerically delicate and at this phenomenological 
stage of our study it does not justify the dedicated effort it requires. For 
certain this critical value indicates the end of the validity of the context of 
our working conditions. This, we expect, indicates a qualitative change on 
the nature of the condensation. 

The region of critical i] is characterised by competing domains of neutral 
and charged scalars, which might present analogies to disoriented chiral con- 
densates jl]. In this region we have estimated the value of r]cr for varying 
self-coupling and present them in Table 1. Although the dependence with A 
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Figure 4: Diagram showing the vacuum selection favouring the charged 
scalars to condense when r/o = in the broken phase. 



is difficult to establish, our results give an indication that as A is increased 
by approximately one order of magnitude, A = 4 — 20, the critical value 
of T] also increases by a similar magnitude, ~ 10^^ — 10~^. This is the 
behaviour we would naively expect from the fact that the scalar self-coupling 
is also a measure of the extent to which the decoupled fields have indirect 
contact withj^he photon bath. Finally, we remark that the results in the 
Table 1 for /5 = 9 might already fall outside the region where a classical 
description based on a Langevin approach is not reliable ^ as the system is 
no longer close to T^,. 

4.3 Interplay between parameters 

It is useful to analyse the interplay between the various parameters in the 
model. To this end, we look at the kinetic energy of each set of fields, 
the coupled and the decoupled. Parallels with the toy model studied in 
the previous section will also be easier to draw. With this study we aim at 
learning how thermalisation is affected by the ratio of the viscosity coefficients 
Tj = tiq/t]^ and the scalar coupling A. 

In Figure 5 the time evolution of the kinetic energies for both types of 
scalar fields are plotted. The parameters r/o = 0.005, r/+ = 1 and [3 = 6.5, 
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Table 1: The table shows the values of (|0i02|) and (|0304|). The three values 
of P correspond respectively to a temperature above, near and below Tc. 
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all in units of the T = vacuum expectation value f in a 100^ box, are the 
same for all the curves, whereas we use three different values for A. 

The equilibrium curve corresponds to the coupled fields which thermalise 
quickly in the time scales displayed in our plot. Of course "quickly" here 
means that the relaxation times for the decoupled fields are clearly longer 
than for the coupled ones. This occurs for all the parameter values shown. 
We verified that the temperature for the equilibrium curve is to a very good 
approximation consistent with the equipartition relation and therefore inde- 
pendent of A. The external dissipation does not give origin to noticeable 
deviations from equipartition, as in the toy model, because of the present 
large number of degree of freedom. 

The most interesting feature of Figure 5 is the A dependence of the asymp- 
totic values for the kinetic energy of the decoupled fields. As in (II Ij) we can 
interpret these quantities as effective equilibration temperatures T^s. We ob- 
serve that the larger A is, the faster the decoupled fields approach a steady 
state. This is to be expected as the decoupled fields interact indirectly with 
the photon bath via the quartic scalar coupling. It explains not only the 
shorter relaxation times for larger values of A but also the higher Tes values 
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Figure 5: Diagram showing the kinetic energy of the scalar fields for a single 
value of ?7 = r/o/^c and three different values of A. For small values of A the 
neutral scalars approaches a steady state at a very slow rate. 



which get closer to the temperature of the bath. 

We verified that the value of the asymptotic kinetic energies is indepen- 
dent of the initial conditions. In Figure 5 the starting kinetic energy of the 
decoupled fields is a third of the equipartition value. It is therefore safe to 
conclude that for weak scalar couplings the equilibrium state should be quite 
distinct from the equipartition thermalised state. 

In Figure 6 we complement the curves shown in Figure 5 by keeping now 
the same scalar coupling for all the curves, here we use A = 8, and vary 
770, or equivalently rj as we again set = 1. We observe that the larger 
the viscosity coefficient r]c the faster the decoupled fields equilibrate. On the 
other hand, as r] increases T^s shifts away from T. The asymptotic steady 
state approaches the equilibrium curve only in the opposite limit, i.e. when 
the external dissipative channel is "switch off". In this case, the kinetic 
energy does eventually reach the value expected by equipartition but at a 
clearly slow rate set by the magnitude of A via the fluctuations mediated by 
the scalar coupling. 
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Figure 6: Diagram showing the kinetic energy of the scalar fields for a single 
value of A and four different values of rj = r]o/rjc. The slowest relaxation time 
for the neutral fields occurs in the r] = limit of no external dissipation. 

4.4 Coupling the decoupled fields 

We end this section with a discussion on the effects of introducing fluctuation 
terms to the evolution equations for the decoupled fields in analogy to what 
we have done in (jl3|) for the toy model. The Langevin equations already 
include indirect interactions of the decoupled fields with the bath. However, 
one may still wonder if results might change when we include a more direct 
source of fiuctuations. The most natural reason for introducing additional 
fiuctuating forces is the interaction between the soft modes described by the 
Langevin equations and the associated scalar hard modes which should be 
present for both types of fields. 

The relative strength of the different interactions is central here. Let g 
denote the gauge coupling. Then our working assumption that the scalar self 
coupling is weak compared to the coupling to the photon bath means g ^ X. 
In a simple perturbative estimate this implies that 77+ ^ t/q. Qualitatively 
this relation is expected to hold at very high temperatures but more reliable 
non perturbative statements require more dedicated simulations. 

In order for the external dissipation to have observable effects it can not 
be negligible with relation to both thermal dissipation coefficients. In our 
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simulations this is trivially the case as 77+ ^ r^g""* = 7]^, but 7]^"^^ > 7]^ = 0. 
The first of these two inequalities justifies neglecting non thermal terms in 
the equations for the coupled scalars. In general, the biasing should occur as 
long as ?7q''* ^ ?7o 7^ 0. Let us see how our analysis supports this view. 

Looking back at equation (fT3j) we see that in the toy model the neutral 
field dissipation must be external for the effective temperatures of the asymp- 
totic states to be different. A similar result applies to the 0(4) model where 
no vacuum selection occurs unless there is an external dissipation. From 
Figure 6 we recognise that it is necessary for the dissipation coefficient 770 to 
have a non thermal external component. Clearly if rjo had a purely thermal 
origin the kinetic energy of the neutral field would always asymptotically 
approach the equilibrium curve. This is expressed in the 77 = curve in 
Figure 6 as thermal dissipative effects are implicitly already present due to 
the scalar self-coupling. Adding an explicit coupling to a heat bath would 
only have the effect of decreasing the relaxation time. 

5 Summary and discussion 

In this paper we discussed the conditions that cause a selective condensation 
of neutral fields in an 0(4) symmetric scalar field sector of a field theory in 
three dimensions when this sector also includes charged fields. In particu- 
lar, we studied a regime in the vicinity of a symmetry breaking transition 
where the interactions with a thermalised photon bath can be simulated by 
a phenomenological Langevin approach. 

We show that under quite general conditions a non thermal dissipation 
in the evolution of the neutral fields effectively reduces the vacuum manifold 
of a system described by an 0(4) scalar model from to S^, above a small 
dissipation threshold. From this analysis we identify the conditions that 
favour the formation and stabilisation of embedded defects as first argued by 
Nagasawa and Brandenberger |T]. 

The vacuum manifold reduction is due to the existence of different asymp- 
totic steady states for the two types of scalar fields considered. This effect 
is caused by an "external" source of dissipation, in the sense that it does 
not arise from fluctuations resulting from interactions with the photon bath. 
The neutral field steady state is characterised by an effective temperature 
Tefj < T, "colder" than the photon bath, and this is what leads to their 
selective condensation. 
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The remarkable feature is that even a small amount of "external" dis- 
sipation can be sufficient to cause qualitatively distinct effects, such as the 
vacuum selection. Small here is in relation to the dominant dissipation terms 
in the charged scalar field sector, which is related to the interaction with 
the photon bath by the fiuctuation-dissipation theorem. This is counterin- 
tuitive. In principle, one would be inclined to neglect the possible effects of 
such a small amount of dissipation. For instance, the asymptotic state of the 
charged field sector is hardly affected by the "external" dissipation. What 
changes things is the existence of a small indirect thermal dissipation. Then 
the external dissipation needs at least to be of the order of the much smaller 
thermal dissipation coming from the indirect coupling of the neutral scalars 
with the heat bath. 

We can naturally generalise the system to a more realistic one. First, 
the external dissipation is considered in both scalar field sectors. However, 
for the charged field sector this leads only to negligible corrections. Second, 
the neutral scalars are coupled directly to the heat bath although with a 
much weaker coupling than the charged scalars so that the resulting thermal 
dissipation for the former does not dominate over the non thermal external 
dissipation. Under these general conditions our results on vacuum selection 
are not qualitatively changed. 

Finally, we note that possible corrections to the scalar fields potential 
coming from the interactions with the gauge bosons should not play an im- 
portant role for our analysis. We know from the work of Nagasawa and 
Brandenberger P that the asymmetry created by the decoupling from the 
neutral fields from the photon bath biases the effective potential in a way 
that stabilises non topological defects when immersed in a photon plasma. 
Moreover, equilibrium thermal corrections to the potential tend to reduce the 
instability of these embedded configurations [21]. Therefore, at a perturba- 
tive level we do not expect corrections to counteract the vacuum selection we 
analyse here. A less investigated difficulty, but potentially an important one, 
is the contribution from very soft photons. Because of infra-red divergences 
reliable corrections similar to those in ^T] ^] are not to our knowledge cur- 
rently available. More dedicated simulations including the full dynamics of 
both the scalars and the gauge bosons are necessary to clarify this problem. 

The type of scenario we describe here should be relevant for applications in 
the early universe where the expansion of the universe provides a non thermal 
source of dissipation. It would also be interesting to investigate if similar 
effects might take place in the quark-gluon plasma where different steady 
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or intermediate states might coexist as in the bottom-up thermahsation in 
heavy-ion colhsion proposed in [20J driven by soft gauge bosons modes. 
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